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Abstract Any constructive continuous function must have a gradually varied approxima- 
tion in compact space. However, the refinement of domain for a — net might be very small. 
Keeping the original discretization (square or triangulation) , can we get some interesting 
properties related to gradual variation? In this note, we try to prove that many harmonic 
functions are gradually varied or near gradually varied; this means that the value of the 
center point differs from that of its neighbor at most by 2. It is obvious that most of the 
gradually varied functions are not harmonic. This note discusses some of the basic harmonic 
. functions in relation to gradually varied functions. 

iO ' 1 Introduction 

In this note, we will discuss some interesting facts about gradually varied functions (GVF) and 
^\ \ harmonic functions. The compatibility between gradually varied functions and harmonic functions 

is important to the applications of gradually varied functions in real world engineering problems. 

Any constructive continuous function must have a gradually varied approximation in compact 
space pp. However, the refinement of domain for a — net might be very small. Keeping the 
original discretization (square or triangulation) , we can obtain some interesting properties related 
to gradual variation. In this note, we try to prove that many harmonic functions are gradually 
varied or near gradually varied, meaning that the value of the center point differs from that of its 
neighbor by at most 2. It is obvious that most of the gradually varied functions are not harmonic. 
This note discusses some of the basic harmonic functions in relation to gradually varied functions. 

Let Ai,A2, A n be rational numbers and A% < A2 < ... < A n . Let D be a graph. / : D — > 
{Ai, A n } is said to be gradually varied if for any adjacent pair p, q in D and f(p) = Ai, then 
f(q) — Ai-i, Ai, or Ai + i. We usually let Ai = i. 

Extending the concept of gradual variation to the function in continuous space: / : D — > R is 
gradually varied if \p — q\ < 1 then |/ g — / p | < 1. Or 

|/,-/p|<|p-?|. (1) 
To some extent, gradual variation is the same as the locally Lipschitz condition. (However, A4 
may be defined differently.) 

On the other hand, a harmonic function satisfies: 

d 2 f d 2 f 

dx 2 ~^ dy 2 ^ ^ 

A main property of the harmonic function is that for a point p, f(p) equals the average value 
of all surrounding points of p. 
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If / is harmonic and p, q are two points such that f(p) < f(q) and s is a path (curve) from p 
to q, then we. We know 

f q -f P = I V/-ds (3) 

If s is a projection of a geodesic curve on /, does the gradient V/ maintain some of its 
properties? For example, is it a constant or does it have any property relating to gradual variation? 
What we would like to prove is that if we define 

fmg{p, q) = max{|V(/)|}on curve s or entire D (4) 

should we have 

Observation A: f m g(p,q) < 2 • \(f q — f p )\/length(s) when / is harmonic? 
Therefore, our purpose is to show that many basic harmonic solutions are at least "near" GVF 
solutions. 



2 Harmonic Functions with Gradual Variation 

Given the value of a set of points in domain D, f : J — > i?,, J C D, for 4-adjacency in 2D (grid 
space), using an interpolating process, we can obtain a GVF solution. pQ We also can solve a linear 
equation using a fast algorithm for a sparse matrix of the harmonic equation based on 

fi,3 = + fi+l,j + fi,j-l + fi.3 + 1) (5) 

or give an initial value for / and then do an iteration. This formula gives a fast solution and 
also gives a definition of discrete harmonic functions [7] . 

How we use the GVF algorithm to guarantee a near harmonic solution is a problem. We can 
use the divide- and- conquer method to have an O(nlogn) algorithm and then iterate it a few times 
to get a harmonic solution. 

Assume 61 and 62 are two points in boundary J. f(bi) < /(^) and s(6i, 62) is a path from b\ 
to b 2 ■ So 

(f(h) ~ f(bt)) 



length(s(bi, 62)) 
is the average slope of the curve. We can define 

slope(h,b 2 ) = max{ ~ { ^\ \s(h, b 2 ) is a path} (7) 

Therefore, there is a s(&i, 62) whose length reaches the minimum. Such a path will be a geodesic 
curve. 

With the consideration of the maximum "slope" , the reason for ObservationA is 

I V/| < (g 2 + |^V /2 <?2 • slope(bl, 62) < 2. (8) 

In general, 

|V/| < + ^ f/ k <?2 ■ slope(bl, 62) < 2. (9) 

where k > 0. Since slope < 1 based on the condition of gradual variation, we want to show that 
the harmonic solution is nearly gradually varied. Note that the gradual variation condition is 
similar to the Lipschitz condition. 

There are two reasons for using "2" in the above formula as the ratio: (1) It is not possible to 
use "1," and (2) anything less than 2 is almost gradual variation. 
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Lemma 2.1 There are simple cases in discrete space that the harmonic solution reaches difference 
1.5. 

Proof Assume that we have five points in grid space in direct adjacency: (i — l,j), (i + 

1, j), - 1), {i,3 + 1) and fi_ ltj = 1, f i+1J = f lJ _ 1 = f iJ+1 = 3 

We want to know what f(i,j) equals. Using the GVF, we get f(i,j) — 2 by Definition 1.1. 
See Fig. 3.1. 
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Figure 1: Discrete harmonic interpolation. 

Using harmonic functions, we will have f{i,j) — 2.5 by Definition 2.1. With the same principle, 
we can let fi-ij — 3 and /i+ij = fi,j-i = fi,j+i = 1- So f(i,j) — 1.5 for the harmonic solution 
and f(i,j) = 2 still for gradually varied. 

When we use the harmonic solution to approximate gradual variation, we need to see if we can 
find the best value when choosing from two possible values. A simple algorithm may be needed 
to make this decision. 

Observation B: There is a GVF that is almost harmonic: \center — averageO f N eighbor\ < 1 
or \center — averageO f N eighbor\ < c, c is a constant. 

The above examples show that a perfect GVF is not possible for a harmonic solution. The 
gradient (maximum directional derivative) less than 2 • f' m , f' m denotes the maximum average 
change (slope) of any path between two points on the boundary possessing the mean of gradual 
variation. 

Every linear function is harmonic. And for quadratic functions, we have 

f(x, y) = ax 2 + by 2 + cxy 
is harmonic if and only if a = —6. However, the following example will not meet the case. 

Example 1 Three vertices of a triangle are pi = (0,0), p2 = (9,0), p3 = (—8,4). The linear 
function f(x, y) = x + 3y. 



(-8.4) 




Figure 2: Example of linear functions. 
This triangle satisfies the gradually varied conditions: 

|/(pl)-/(p2)|=9<|pl-p2|=9 

|/(p2)-/(p3)| = |9-4|<|pl-p2| 
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|/(pl)-/(p3)| = |0-4| < V8 2 + 42 

If we consider a point p — (x, y) on the line < P2,P3 > when x = and y = 36/17 = 2.4, then 
f(x, y) > 7. This point and p\ do not maintain the condition of gradual variation. \f(p) — f(pi)\ > 
7 > \p-pi\. 

This example seems to break the observation we have made. However, let us revisit the function 
f(x, y) = x + 3y and let z = f(x, y). We have z — x — 3y = 0. We can have y = ^z — |x represent 
the triangle and associated function. In general, for a linear function in 3D 

ax + by + cz + d = 

We can always find a coefficient that has the maximum absolute value. We will have the 
equivalent equation that has 

AX + BY + D = Z (10) 
where \A\ and \B\ < 1. This property is often used in computer graphics. 

Lemma 2.2 The Piece-wise linear function preserves the property of gradual variation. 

Proof: We first want to discuss the case of a single triangle where any piecewise linear function 
is a harmonic function. In this case we can write the function like this 

f(x,y) =ax + by + c, |a|,|6|<l 

^ = f x = a , |£ = f v = b. The gradient is a constant \J a 2 + b 2 . There is a horizontal and 
vertical line that goes through boundary points. The maximum average rate of change r (average 
slope on the path between two points on the boundary) is greater than or equal to max a, b. 

Since \a\, \b\ < 1 ; r < Va 2 + b 2 < \/2maxa,6 < s/2. So r < 2. 

If this piecewise linear function is on a polygon (2D), it will still have this property. 

The problem is that in this proof, we have not used the conditions of gradual variation directly. 
The conditions are 

|/(pi) - f(P2)\ = \a(xl - x2) + b(yl - y2)\ < \ Pl - p 2 \ = y/(xl - x2) 2 + (yl - y2) 2 
l/(pi) - /(Ps)l = Hxl - x3) + b(yl — 2/3) j < \ Pl - p 3 \ - ^ (xl - x3) 2 + (yl - y3) 2 

|/(P2) - /(Ps)| - Hx2 - x3) + b(y2 - y3)\ < \p 3 - p 2 \ = \[(x2 - x3) 2 + (y2 - y3) 2 

We have used the GVF general property and the triangle constraint. The next section will 
discuss a more general case. 



3 Gradually Varied Semi-Preserving 

In this section, we extend the content of above sections using more rigorous mathematical def- 
initions. Harmonic functions can be characterized by the mean value theorem. Here we are 
interested in harmonic functions that are gradually varied. More specifically, a function is said to 
be gradually varied semi-preserving if 
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i^-, i \ u (p) — u (q)\ /-, , x 

max Vu <c- max J— ^ ^ (11 

where Vu is the gradient of u, D is a domain with the boundary 9D, and c is a constant. 

The above formula poses a property of computational importance. We can show that lin- 
ear functions and quadratic hyperbolic functions satisfy the condition of gradually varied semi- 
preserving. 

If u is linear we can assume that u = ax + by + c and if u is quadratic hyperbolic we can let 
u = a(x 2 — y 2 ). We do not restrict the value of a, b, c here. 

Lemma 3.1 Proposition 2 If u is linear or quadratic hyperbolic, then 

max|V U |<y(2). max = (12) 

B v v / p qedB \p-q\ 

where B is any ball. 

Proof Let u be a linear function u — ax + by + c then 

|V«| - V(a 2 + b 2 ) < v^2) max{|o|, |6|} (13) 

On the other hand, if we choose p — (— r, 0), q = (r, 0) on dB, 
where r is the radius of the ball B. Then 

\u(p)-u(q)\ _ \-ar-ar\_ = 



\p — q\ 2r 
Choosing another pair of p and q on dB, 

p= (0,r),g=(0, 

we have 



Hp)-u(q)\__ lbl 



\p-q\ 

Combining (13), (14) and (15) we conclude (12) when u is linear. 

Now, consider u as a quadratic hyperbolic function: u = a(x 2 — y 2 ). Then, 

|Vu| = 2\a\ v r (x 2 +y 2 ) < 2\a\r (16) 
On the other hand, if we choose p and q on dB, 

p=(0,r),g=(r,0) 

Then 

\u(p)-u(q)\ 



\p-q\ vT 2r2 ) 

Combining (16) and (17) we have 



2\a\r (17) 



vu| < V( 2 )Mp)^M < 7(2) ~~ w*) - u ^ 



max 



\p — q\ x,y£dB \x — y 

and (12) follows. 
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4 Discussion 



Recent studies show an increased interest in connecting discrete mathematics with continuous 
mathematics, especially in geometric problems. For instance, the variational principle has been 
used for triangulated surfaces in discrete differential geometry, see [5] . This note presented an idea 
of combining a type of discrete function: the gradually varied function and a type of continuous 
function: the harmonic function in a relatively deep way in terms of continuous mathematics. 
The harmonic function is a weak solution to the Dirichet problem which is about how to find a 
surface when the boundary curve is given. The gradually varied function was proposed to solve 
a filling problem in computer vision. We are hesitant to use the method of the Dirichlet prob- 
lem for the discrete filling problem since we do not know the exact formula (function) on the 
boundary, even though we know the sample points. This problem is also related to the Whitney's 
problem 5 6]. Some ideas have been presented by the first author in the Workshop on Whit- 
ney's problem in 2009 organized by C. Fefferman and N. Zobin at College of William and Mary 
(http://nxzobi.people.wm.edu/whitney/whitney.htm). During the Workshop, P. Shvartsman pre- 
sented an idea of using geodesic curves in Sobolev space (See related paper [9]). Our idea about 
using the geodesic curve presented in this note was independently obtained. 
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